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ABSTRACT 

Context. The blazar 1 ES 1218+30.4 has been previously detected by the VERITAS and MAGIC telescopes in the very high energies. 

The new detection of VERITAS from December 2008 to April 2009 proves that 1 ES 1218+30.4 is not static, but shows short-time 

variability. 

Aims. We show that the time variability may be explained in the context of a self-consistent synchrotron-self Compton model, while 

the long time observation do not necessarily require a time-resolved treatment. 

Methods. The kinetic equations for electrons and photons in a plasma blob are solved numerically including Fermi acceleration for 

electrons as well as synchrotron radiation and Compton scattering. 

Results. The light curve observed by VERITAS can be reproduced in our model by assuming a changing level of electron injection 

compared to the constant state of 1 ES 1218+30.4 . The multiwavelength behaviour during an outburst becomes comprehensible by 

the model. 

Conclusions. The long time measurements of VERITAS are still explainable via a constant emission in the SSC context, but the short 

outbursts each require a time-resolved treatment. 
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1. Introduction 

Blazars are a special class of active galactic nuclei (AGN) 
exhibiting a spectral energy distribution (SED) that is strongly 
dominated by nonthermal emission across a wide range of 
wavelengths, from radio waves to gamma rays, and rapid, 
large-amplitude variability. The source of this emission is 
presumably the relativistic jet emitted at a narrow angle to the 
line of sight to the observer 

In high-peaked BL Lac objects (HBLs) the SED shows a double 
hump structure as the most notable feature with the first hump in 
the UV- to X-ray regime and the second hump in the gamma-ray 
regime. Indeed, a substantial fraction of the known nearby HBLs 
have already been discovered with Cherenkov telescopes like 
H.E.S.S., MAGIC or VERITAS. The origin of the first hump 
is mostly undisputed: nonthermal, relativistic electrons in the 
jet are emitting synchrotron radiation. The origin of the second 
hump is still controversially debated. Up to now two kinds of 
models are discussed: leptonic (e.g. Maraschi et al. 1992 1 and 
hadronic (e.g. Mannheim 1993) ones, which are mostly applied 
for other subclasses of blazars. 

Another important feature of AGNs in general and HBLs in 
particular is their strong variability. The dynamical timescale 
may range from minutes to years. This requires complex 
models, which obviously have to include time dependence, but 
this gives us also the chance to understand the mechanisms that 
drive AGNs. We will apply a self-consistent leptonic model to 
new data observed for the source 1 ES 1218+30.4, because 
those are the ones favoured for HBLs. 



The source HBL 1 ES 1218+30.4 has been discovered as a 
candidate BL Lac object on the basis of its X-ray emission 
and has been identified with the X-ray source 2A 1219+30.5 
(Wils on et al. 1979; Led den et al.|[T981 ). For the first time, 
1 ES 1218+30.4 has been observed at VHE energies using the 



MAGIC telescope in January 2005 (lAIbert et al.||2006| ) and 
later from VERITAS ( |Acciari et aL]|2009| l. Co verage of the 
optical/X-ray regime is provided by BeppoSAX (Donato et al. 
2005) and SWIFT (|T ramacere eT aT 2007), unfortunately the 
data are not always simultaneous. During the observations 
from December 2008 to April 2009 VERITA S also observed 
1 ES 1218+30.4 showing a time-variability ( |the VERITAS] 
collabor ation et aLl|2010| ). The observations from the MAGIC 
telescope have previously been modelled by Riiger et al. (2010 1. 
[the VERITAS collaboration et al] ( |201()| ) claimfeat their 
new observations exhibiting variability challenge the previous 
models. We will show that a timedependent model using a 
self-consistent treatment of electron acceleration is able to 
model the new VERITAS data. 



We present the kinetic equation, which we solve numeri- 
cally, describing the synchrotron-self Compton emission (Sect. 
|2]). In Sect. [3] we apply our code to 1 ES 1218+30.4, taking 
the VERITAS data into account and give a set of physical 
parameters for the most acceptable fit. Finally, we discuss 
our results in the light of particle acceleration theory and the 
multiwavelength features. 
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2. Model 

Here we will give a brief description of the model used, for 

a complete ov erview see ( [Weidinger et al.|[2010[ [Weidinger &| 

Spanier|2010| l. 

We start with the relativistic Vlasov equation (see e.g. 



Schlickeiser 



2002 1 in the one dimensional diffusion approxima- 



tion (e.g. Schlickeiser|1984[ l, here the relativistic approximation 
p X ymc is used. This kinetic equation will then be solved 
time-dependently in two spatially different zones, the smaller 
acceleration zone and the radiation zone, which are assumed to 
be spherical and homogeneous. Both contain isotropically dis- 
tributed electrons and a randomly oriented magnetic field as 
common for these models. All calculations are made in the rest 
frame of the blob. 

Electrons entering the acceleration zone (radius 7?acc) from the 
upstream of the jet are continuously accelerated through diffu- 
sive shock acceleration. This extends the model of IKirk et al.l 
( |199 8 ) with a stochastic part. The energy gain due to the accel- 
eration is balanced by radiative (synchrotron) and escape losses, 
the latter scaling with fesc = '7^acc/c with rj - 10 as an empiri- 
cal factor reflecting the difiiisive nature of particle loss. Escaping 
electrons completely enter the radiation zone (radius /?iad) down- 
stream of the acceleration zone. 

Here the electrons are suff'ering synchrotron losses as in the 
acceleration zone and also inverse-Compton losses, but they 
do not undergo acceleration. Pair production and other contri- 
butions do not alter the SED in typical SSC conditions and 
are neglected ( [Chiang & B ottcher 2002). The SED in the ob- 
server's frame is calculated by boosting the selfconsistently cal- 
culated photons towards the observer's frame and correcting for 
the redshift z: ly^^^ = 5^/iVob.s/(47r)A^ph with Vobs = <5/(l + z)v. 
The acceleration zone will have no contribution to /obs directly, 
due to the R} dependence of the observed flux at a distance r 
(Fy^^^ir) - jiIy^^^R^_^^r^^) and the small size of the acceleration 
zone. The kinetic equation in the acceleration zone is 



^=|;[(^./-Cr)-«.(r,0]. 



d_ 
dy 



[(a + 2)facc] r 



dy 



+ Qq(7 - ro) - tjMr' t) . 



(1) 



The injected electrons at yo, as the blob propagates through the 
jet, are considered via Qmjiy,!) :- QoS(y-yo)- The synchrotron 
losses are calculated using Eq. (|2]l. 



Ps(y) 



1 o-tB^ 2 



67T mc 



7'-/3sr 



(2) 



with the Thomson cross-section ctj. The characteristic timescale 
for the acceleration face - (^i/(4^||) + 2v^/(9A'||)j of the sys- 
tem is found by comparing Eq. ([TJ with [Schhckeiser[ ( 1984| l with 
the parallel spatial diffusion coefficient K\\ not depending on y 
when using the hard sphere approximation. The characteristic 
timescale has an additional factor (oc y^) arising from the 
Fermi-II processes compared to shock acceleration by itself. 
The stochastic part of the acceleration also gives rise to the 
second row in Eq. ([T}, while the first row mainly depends 
on Fermi-I processes. This dependence of face is important 
for the interpretation of the resulting electron spectra, e.g. of 
their slopes (depending on facc/fesc) or the maximum energies 
(depending on 1 /(facc/?i)), see Weidinger et al. ( 2010J ) for details. 
For modelling SEDs and lightcurves it is primary important to 



ensure sensible values for face- Unlike in [Drury et al.[ ( [1999| l, 
the energy-dependence of the escape losses is also neglected 
because we do not expect a pileup as suggested in S chlickeiser] 
( 1984 1 at typical SSC conditions, v^, va are the shock and Alfven 
speed respectively. Hence a in Eq. ([!]) measures the efficiency of 
the shock acceleration compared to stochastic processes. Setting 
va - 0, i.e. a — > oo, will result in a shock-only model like |Kirk[ 
^taL] ( fT998l l. 

This model takes account of a much more confined shock 
region. Fermi-I acceleration will probably not occur over the 
whole blob when considering a real blazar but rather at a small 
region in the blob's front. Neglecting acceleration simplifies the 
kinetic equation in the radiation zone to 



dNeiyJ) d 



dt 



dy 
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Pic accounts for the inverse-Compton losses of the electrons 
additionally occurring (beside the synchrotron losses) (e.g. 
|Schhckeiser|2002| i: 



Picij) ^m c 



h I daa I 

Jo Jo 



daxN^kiax) 



dN(y, ai ) 
dtda 



(4) 



The photon energies are rewritten in terms of the electron rest 
mass, hv = amc^ for the scattered photons and hv - a\mc^ 
for the target photons respectively. Equation (|4]l is solved nu- 
merically using the full Klein-Nishina cross-section for a single 
electron scattering off a photon field (see e.g. Jones 1968 ). Here 
a,„ax accounts for the kinematic restrictions on IC scattering. In 
analogy to the acceleration zone the catastrophic losses are con- 
sidered via fescrad - ^Rrdd/c with ?7 = 10. fesciad IS the respond- 
ing timescale of the electron system, which is proportional to the 



variability timescale in the observer's frame (see e.g. Kerrick et 
[aL|r995i ): 



' var ^ 



(5) 



To determine the time-dependent model SED of blazars the par- 
tial differential equation for the differential photon number den- 
sity has to be solved time-dependently, which can be done nu- 
merically. The PDF (|6| can be obtained from the radiative trans- 
fer equation making use of the isotropy of the blob 



'9A^ph(v,f) ^ „ , , „ ^Ph(^'0 



dt 



(6) 



'ph.esc 



where Rs and Re are the production rates for synchrotron photon 
and the inverse-Compton respectively. Rs is calculated using the 
well known Melrose approximation and the inverse-Compton 
production rate Re is treated in the most exact way, i.e. using 
the full Klein-Nishina cross section, see Weidinger et al. (2010}. 
Below a critical energy the obtained spectrum is self-absorbed 
due to synchrotron self-absorption, which is described by ay 
( [Weidinger et al.|2010[[Ruger et al.|20101 l. The photon-loss rate 
is set to be the light-crossing time. 



3. Results 

Using the parameters summed up in Table [T] we were able to 
fit the emission of 1 ES 1218+30.4 as a steady state with our 
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SSC model, see Fig. 1. We used all the archival data from 
BeppoSAX, SWIFT in the X-ray band and the MAGIC 2006, 
VERITAS 2009 as well as the new released VERITAS 2010 



data in the VHE to model the SEP of 1 ES I2I8+30.4 (|Donato 
eraL||2005| [T ramacere et al.ll2007', l AJbert et al.||2006[ 



etal.|2009 



Acciari 



the VERITAS collaborati on et al.|2010| l. Th e derived 



SED is absorbed in the VHE using the EBL model of'Primack' 
etaL] (p005 ) for the corresponding redshift of 1 ES 1218+30.4 . 
The parameters of our SSC model are well winthin the standard 
SSC parameter region with an equipartition parameter of 0.02. 
Even though PIC and MHD simulations suggest a higher 
magnetic field compared to particle energy (in the range of 0. 1), 
this is a common assumption in SSC models, but has to be kept 
in mind with regard e.g. the stability of the blob. If one wishes 
to enforce higher equipartition parameters one could to use the 
model of Schlickeiser & Lerche (2007). In order to allow strong 
shocks to form va < vs must be fulfilled, which is the case for 
fl = 10. 

Due to relatively small deviation (within the error margins) 
between the MAGIC 2005/VERITAS 2008 and the averaged 
VHE data from the VERITAS 2009 campaign we find a steady 
state the most plausible way to model the emission, i.e. the small 



fluctuations (see the overall lightcurve in the VERITAS collab- 
oration et al. (2010)l) are not contributing significantly to the 



averaged observed SEDs. In the Fermi LAT energy regime our 
model yields a photon index of apei = - 1 .69, whichagrees weU 
with the Fermi measurement of - 1 .63 + 0. 12 ( |Abdo et al.|2009| l. 
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Fig. 1. Model SED of 1 ES 1218+30.4 (black soHd line) as de- 
rived using the described model (see Sect. l2]i and the parame- 
ters shown in Table [T] The VHE parts of the model SEDs have 
been absorbed using the EBL model of |Primack et al. (2005|l. 



The BeppoSAX data are from Donato et al. 



from jTramacere et al.l(|2007ll, MAG IC from'Al bert et al.| ( |2006l l 



(|2005|l, SWIFT 



VERITAS 2009 from |Acciari et al.|p0 09) and the blue dots are 
the new VERITAS 2010 data from [fee VERITAS collaborationi 



et al. (2010). The dashed red curve shows the time integrated 
SED over the strong outburst shown in Fig. 2, as measured by 
VERITAS in 2009. 



Table 1. Model parameters for the low-state SED, basis to model 
the outburst by varying Qo. 



eo(cm-^) 


fi(G) 


^acc(cm) 


^iad(cm) 


'acc / ^esc 


a 


6 


6.25 • 10" 


0.12 


6.0 ■ 10'* 


3.0-10'= 


1.11 


10 


44 



The lightcurve of [the VERITAS collaboration et al.| (J20T0]l 



shows a relatively strong outburst at =s MJD5486I. Starting with 
the steady state emission (solid line. Fig. 1; parameters: Table[T]l 
we injected more electrons Qo into the emission region at low 
yo ~ 3. As the blob evolves in time the emission in the model 
at higher energies rises and drops off again when the injected 
electrons finally relax to the initial Qq. This process can be 
explained as density fluctuations along the jet axis and finally 
fits the flare. 

We found that nearly doubling the injected electron number 
density in a Qc(i) = 1 + bitlt^^^^^f' way with a timescale 
fe.var ~ 1 .5 days (as measured in the observer's frame) and then 
decreasing them to the initial Qq in an almost linear way on the 
same timescale fits the strong outburst of 1 ES 1218+30.4 . The 
coiTesponding lightcurve of the model as well as the observed 
one are summarized in Fig. 2 3. Figure 3 shows a more detailed 
view of the lightcurve in the VHE (above 200 GeV) as well as 
the corresponding lightcurves in the X-Ray (between 1.2 keV 
and 1 1 keV) regime of BeppoSAX/SWIFT and the lower tail of 
the Fermi LAT energy range (between 0.2 GeV and 22 GeV) as 
predicted by our model. The latter two have been scaled down to 
the flux level of the VERITAS measurement, see Fig. 3, because 
the real fluxes are higher than the VHE flux. The model predicts 
the peak of the lightcurve in the Fermi regime to be 1 .26 hours 
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Fig. 2. Lightcurve of the photon flux above 200 GeV as mea- 
sured by|flie VERITAS collaboration et al."'('20T0) (inset of their 



figure) in January 2009 to February 2009 and our model (red 
solid line). The outburst was modelled by injecting more elec- 
trons into the blob by varying Qq at a timescale of ^ 1.5 days 
(see text for details). 



ahead of the VHE one, where the X-ray regime is delayed by 
0.97 hours for 1 ES 1218+30.4 . The delay of the X-ray band 
can be used to verify the model when multiwavelength data of 
the flaring behaviour of 1 ES 1218+30.4 is available, while the 
derivation of the 0.2 GeV to 22 GeV lightcurve is beyond the 
resolution of Fermi for this source. 

Additionally we plotted the time averaged SED (over the 
outburst from MJD54860 until MJD54864) into the SED of 
1 ES 1218+30.4, Fig. 1 (dashed red line). As one can see 
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our SSC model, I(E > 200 GeV) 
VERITAS 2010. 1(E > 200 GeV) 
model: 1(0.2 GeV < E < 22 GeV)/(1.21xl0"') 
model: 1( 1.2 keV < E < 11 keV)/(2.74xlo"l 




6 7 
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Fig. 3. Detailed view of the high outburst shown already in Fig. 
2 as well as the behaviour of 1 ES 1218+30.4 in the lower Fermi 
LAT energy band and the synchrotron regime, measurable by 
BeppoSAX/SWIFT during such a flai-e. 



only when separately considering the strongest outburst of 
1 ES 1218+30.4 within the VERITAS campaign in 2009 the 
aberration from a presumed steady state is significant. In 
contrast an average over the whole observation of |the VERITAS 
Icollaboration et aL] ( |2010| l, which is low-state most of the time, 
will re sult in a steady state as shown here or in .Riiger et al. 



( |2010| l. For the IC photon index a (vFy oc v«+2) above 200 GeV 



we get cKvHE - -3.53 for the low-state (solid curve in Fig. 1), 
which slightly softens to q-vhe = -3.56 when considering the 
high-state as the time-average over the single outburst shown 
in the lightcurve. Fig. 2 in the VERITAS range. Note that 
the photon index and its behaviour during an outburst in this 
energy range is very sensitive to the EBL absorption and thus 
to the EBL model used and shows a strong dependency on 
the considered energy range. With our model we are able to 
compute the spectral behaviour in the X-Ray energy range of the 
BeppoSAX/SWIFT satellites (i.e. 1.2 keV < E < 11 keV). The 
model predicts the source to be spectrally steady in this regime 
with a photon index a- 



xray 



-2.68 for outbursts on timescales of 
days. Considering shorter averaging timescales of the outburst 
of 1 ES 1218+30.4 , e.g. the first or last two hours, two hours 
around the peak in the lightcurve, the maximum derivation from 
Cxray (^vhe) predicted by the model is ±0.05 (-0.07), which 
could not be measured with current experiments and thus is 
considered as spectrally steady in this case. 



4. Discussion 

Our results clearly show that the latest observations from the 
VERITAS telescope for 1 ES 1218+30.4 still agree with a con- 
stant (steady state) emission from a SSC model when averaged 
over a long observation period. This is due to the relatively mod- 
erate variability of 1 ES 1218+30.4 compared to the observation 
time. 

The variability may be well explained in the context of the self- 
consistent treatment of acceleration of electrons in the jet. We 
are aware that an outburst of the timescale of roughly five days 
as measured from 1 ES 1218+30.4 does not necessarily require a 
shock in jet model, which scales down to a few minutes depend- 
ing on the SSC parameters ( Weidinger & Spanier 20 10), but may 
also be explained as e.g. different accretion states. Nevertheless 
the fundamental statement remains the same: long time obser- 



vation of slightly variable blazars will result in a steady state 
emission, while an average over a single outburst will, of course, 
result in a significantly diff'erent SED for the source. We are 
not yet able to rule out different emission models or even com- 
plex geometries of the emitting region. But we are able to model 
the influence of short outbursts of a source on the SED and the 
lightcurves in the different energy bands selfconsistently. 
The VERITAS collaboration only shows an integrated spec- 
trum for 1 ES 1218+30.4 , which is due to the low flux of the 
source and the photon index behaviour of the combined high- 
states. This integrated spectrum does not show strong variations 
with regard to the known low-state observed by MAGIC. Our 
model now predicts a clear change in the spectrum, which is in- 
dicated by the dashed line in Fig. 1, which shows the average 
over one outburst with a slight, currently not detectable spec- 
tral softening in the VHE range, while the synchrotron peak 
in the BeppoSAX/SWIFT regime remains spectrally constant. 
This situation changes for shorter and/or stronger outbursts of 
an overall timescale of hours, which will result in measurable 
spectral evolutions in all energy regimes when considered with 
the presented model. Furthermore the time-resolved SEDs dur- 
ing a flare are comprehensible with our model. Hence with better 
time-resolved spectra or/and better multiwavelength coverage it 
should be possible to prove this model, and if the model is in- 
deed applicable it will be a good tool to investigate the whole 
SED during an outburst without having all energy regimes ob- 
servationally covered. 
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